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_ 1 Introduction 

o , 

G^ , We consider evolutionary equation 



O 



X 



Ut + UU^- K {u'^Ux)^ = {u-Uij Lp{u), (1) 

where k, Ui are positive constants, ^{u) is a smooth function nullifying at the origin 
and maintaining a constant sign within the set (0, kUi) for some k> 1. Equation ([1]) 
referred to as convection-reaction-diffusion equation is used for simulating transport 
phenomena in active media. Due to its practical applicability and a number of unusual 
features, equation ([1]) was intensely studied in recent decades [Il[2l[3llll[5l[6l[7l[8]. 

Present investigations are mainly devoted to qualitative and numerical study of the 
family of traveling wave (TW) solution to equation ([T]). Our aim is to show that under 
certain conditions the set of TW solutions contains periodic regimes, solitary waves and 
compactons with special emphasis put on the existence of compactons. Solitary wave 
solutions (or solitons) are exponentially localized wave packs moving with constant 
velocity without change of their shape. They are mostly associated with the famous 
Korteveg de Vries (KdV) equation and members of KdV hierarchy having the form 

K{m) = Ut + v^Ux + Uxxx = 0, m > 1. (2) 
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For m = 1 solitary wave solution to KdV is as follows [9]: 



u{t, x) = U{0 = U{x -Vt) 



cosh (^^e 

where V stands for the wave pack velocity. 

In 1993 Philip Rosenau and John Hyman put forward the following generalization 
to KdV hierarchy [10]: 



K{m, n)=ut+ {un. + K)... = 0, m,n>2. (3) 

Modification in the higher- derivative term causes that members of K{m, n) hierarchy 
possess solutions with compact support. For m = n = 2 compactly supported solution 
(or compacton) has the form 



u{t, x) = f/(0 = U{x-Vt) 



iV cos2 [2^] if 

otherwise. 



x~Vt\ 



Both solitons and compactons are the TW solutions, depending in fact on a single 
variable = x — V t, are described by ODEs appearing when the TW ansatz u{t, x) = 
U{^) is inserted into the source equations (for details see e.g. [H]). Therefore it 
is possibile to give a clear geometric interpretation of solitons and compactons. To 
both of them correspond homoclinic trajectories in the phase space of the factorized 
equations. The main difference between the solitary wave and compacton is that the 
first one is nonzero for any ^ while the last one is nonzero only upon some compact set. 
This is because the soliton corresponds to the homochnic loop usually bi-asymptotic 
to a simple saddle. Hence the "time" needed to penetrate the close loop is infinite. 
Compacton corresponds to the homochnic loop bi-asymptotic to a topological saddle 
lying on a singular line. A consequence of this is that the corresponding vector field 
does not tend to zero as the homoclinic trajectory approaches the stationary point, 
hence the "time" required to penetrate the trajectory remains finite. Let us note 
that the compacton in fact is a conjunction of the nonzero part corresponding to the 
homoclinic loop and the trivial constant solution represented by the saddle point [TT] . 

For the sake of brevity we maintain the notation traditionally used in more specific 
sense. Thus, soliton is usually associated with the localized invariant solution to any 
completely integrable PDE, possessing a number of unusual features [9]. Some of these 
features are inherited by compactons as well |10l [121 [T3| . Here we maintain this notion 
to those solutions to ([1]), which manifest similar geometric features as "true" wave 
patterns known under these names, not pretending that they inherit all features of 
their famous precursors. 

The structure of the article is following. In section 2 we present the dynamical 
system describing the TW solutions to ([1]) and make its qualitative analysis, obtaining 
conditions that contribute to the homoclinic loop appearance. At the end of this 
section we formulate the conditions enabling to make a distinction between compactons 
and solitons. Conditions formulated in section 2 contribute but not guarantee the 
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homoclinic loop appearance. The following section contains the results of numerical 
simulations revealing that expected scenarios of solitons' and compactons' appearance 
really take place. In section 4 we perform a brief discussion of the results obtained and 
outline directions of further investigations. 

2 Factorized system and its qualitative analysis 
2.1 Statement of the problem 

We are going to analyze the set of TW solutions to ([T]), described by the following 
formula: 



Inserting ansatz (jl]) into the GBE one can obtain, after some manipulation, the fol- 
lowing dynamical system: 



where A(f/) = kU". We are going to formulate the conditions contributing to soliton- 
like and compacton-like solutions to equation ([1]) by analyzing the factorized system 
([5]). Let us note that in the case of both KdV and K{m, n) hierarchies' members 
the procedure of factorization lead to the Hamiltonian dynamical systems for which 
distinguishing the presence of homoclinic loops is more or less trivial [H]. The situation 
with system ([5]) is not so simple since it is not Hamiltonian. And we cannot expect 
that the homoclinic loops will form a one-parametric family as this is the case with the 
Hamiltonian systems. Contrary, the homoclinic loop in our case can appear as a result 
of one or several successive bifurcations taking place at specific values of the parameters. 
Let us stress that existence of the homoclinic loop corresponding to compacton is 
possible due to the presence of the factor A(f/) at the LHS of system Since the 
function A{U) contains the origin, which is the stationary point of system ([5]), then 
the natural way of proceedings is following. We state the condition that guarantee the 
appearance of the stable limit cycle in proximity of the stationary point (f/i, 0), and 
ensure that simultaneously the stationary point (0, 0) is a saddle and remain so as the 
parameter of bifurcation is changed. By proper change of the bifurcation parameter we 
can cause the growth of size of the limit cycle and in presence of near-by saddle it finally 
could give way to the homoclinic bifurcation. Since the latter bifurcation is nonlocal, 
we cannot trace its appearance using the methods of local nonlinear analysis, as this 
is the case with Andronov-Hopf bifurcation. Therefore on the final step we are forced 
to resort to numerical simulations. And of course answer the question on whether the 
homoclinic loop corresponds to compacton-like solution or not needs special treatment 
based upon the asymptotic analysis. This and other issues are analyzed in the following 
subsections. 




(4) 



A{U) U = A{U) W, 
A{U) W = ip{U) {Ui -U)- KnU""-^ W^ + {U - V) W, 



(5) 
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2.2 Andronov-Hopf bifurcation in system (|5]) 

To formulate the conditions which guarantee the stable limit cycle appearance in vicin- 
ity of the stationary point (f/i, 0), let us consider the Jacobi matrix 

A(f/i) 
-v(Ui) Ui-V. 

corresponding to this point. In order that {Ui, 0) be a center, the eigenvalues of Ji 
should be pure imaginary and this is so when the following conditions are fulfilled: 

Trace J^ = Ui-V = 0, (6) 
Det Jl = A(t/i)^(f/i) > (7) 

The first condition immediately gives us the critical value of the wave pack velocity 
Vcr = Ui. The second one is equivalent to the statement that v^(f/i) is positive. 

The next thing we are going to do is a study of the stability of limit cycle appearing. 
As is well known [El [15], this is the real part of the first Floquet index 3ftCi that 
determines whether the cycle is stable or not. 

For A{Ui) = kUi > condition 3ftCi < assures that the limit cycle, appearing in 
system when V < Vcr, will be stable. 

To obtain the expression for 3ftCi, the standard formula contained e.g. in [15] could 
be applied providing that system in vicinity of the center is presented in the form 

il \ ^ / -n\ ( Zi\ ( F{zu Z2) 
Z2 ) \^ )'\Z2 G{ZuZ2) 

where VL = y/fj, ■ u, fj, = A{Ui), u = f{Ui), F{zi, Z2) and G{zi, Z2) contain all nonlinear 
terms. In this (canonical) representation the real part of the Floquet index 3ft Ci is 
expressed by the following formula p^ : 

16 9ft Ci = Fill + -^122 + 6*112 + 6*222 + ^ {-^12 (-^11 + -^22) " 

—6*12 (Gil + G22) — Fii Gil + -^22 G22} • (9) 

Here Fjjfc, Fij stand for the coefficients of the function's F{zi,Z2) monomials ZiZjZk-, 
Zi Zj correspondingly. Similarly, indices Gijk-, Gij correspond to the function's G{zi, Z2) 
monomials' coefficients. 

Assuming that relations ([6])-([7]) are satisfied, let us rewrite system ([5]) in coordinates 
y^ = U~Ui, y2 = W: 

u \ f y^ \ . f $i(yi, y2) 



yi 



^^^^d^ \y2 ) J \ y2 J ' \ ^2iyi, ^2) 



where $i(yi, 1/2), $2(2/1, 2/2) are nonlinear terms which are as follows: 

1 



^i{yi, 2/2) = f^nU"^ 'yiy2{Ui + 



n 



$2(2/1, 2/2) 



{KurV2 + 

+2/1 



n[n 



2 

m-2„ 



yi]+o{M') 



l)«:f/rV2/2 + 



(10) 

(11) 

(12) 



viUi) y^ + -(p{U{) yl 



2/2 



o (i2/.r) 



A passage to canonical variables {zi, Z2) can be attained by means of transformation 



Z2 



-v^ 
^ 



yi 
y2 



This gives us the canonical system ([8]) with 



G{zi,Z2] 



F(ZI,Z2] 



2 



n 



-Zl Z2 



n — 1 



-Zl 



+ Oi\z.\'), (13) 



+ ^^^^zl [in -l)zi- UiV^ + Oi\z,\% 



2un 



n 



Analysis of the formulae ([9]), f|T3|l shows that the function F(zi,Z2) does not con- 
tribute to the real part of the Floquet index, which in our case is expressed as follows: 



16Q 



G12 (Gil + 6*22) — — 



len^ifiUi 



{A{Ui) ^(Ui) + Kn^{Ui)Ur'} ■ 



So the following statement holds. 
Theorem 1. If the inequality 



Ui ifiUi) + rnpiUi) > 



(14) 



is fulfilled then in some vicinity of the critical value of the wave pack velocity Vcr = Ui 
a stable limit cycle exists. 



2.3 Study of the stationary point (0, 0) 

From now on we shall investigate system with ^{U) = U"^: 

A{U) if = A{U) W, 
A{U) W = ip{U) {Ui-U)-Kn U""-^ W^ + {U- V) W, 



(15) 



Our next step is to state the conditions assuring that the stationary point (0, 0) is 
a topological saddle or at least contains a saddle sector in the right half-plane. The 
standard theory ^16j can be applied for this purpose. Our system can be written down 
in the form 



d 

df 



U 

w 





Ui 0(0) 





-V 



u 
w 



+ nonl. terms, 



(16) 



where = k U"' ^ and since the trace of its linearization matrix is nonzero then the 
analysis prescribed in Chapter IX of [16] for this type of a complex stationary points 
is the following. 
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1. Find the change of variables {U, W, T) i— {x, y, r) enabhng to write down system 
( IT6|) in the standard form 

-r = P2 (a;, y) , 

= y + Q2ix,y), 

a T 

where P2 (x, y), Q2 (x, y) are polynomials of degree 2 or higher. 

2. Solve the equation y + Q2 {x, y) = with respect to y, presenting result in the 
form of the asymptotic decomposition y = if{x) = ai x^^ + 0:2 + .... 

3. Find the asymptotic decomposition 

P2(x, ^(x)) = A^x'" + .... 

4. Depending on the values of m and the sign of A^, select the type of the complex 
stationary point using the theorem 65 from [16]. 

5. Return to the original variables ([/, W, T) and analyze on whether the geometry 
of the problem allows for the homoclinic bifurcation appearance. 

So let us present results obtained for system (fT5|) . The case m > 1 is most simple for 
analyzing since the canonical system is obtained by the formal change (t/, W) ^ (x, y) 
and passage to new independent variable r = —VT. As a result we get the following 
system: 

dt tjc - ^ , ^ 

— = --x^y=P2{x,y), (17) 
a r V 

4^ = y - [xy - nnx""'^ y"^ - x"" (x-f/i)} = y + Q2{x, y). 

CLT V 



Presenting y in the form of series y = aiXi+... and solving the equation y + Q2{x, y) 

Hi 

V 



0, we obtain y = Lp(x) = ^x^ + .... Inserting function (p{x) into the RHS of the first 



equation we get 

P2{X, ip{x)) = A„+„X"+™ + ... = -K f/i X^^"' + .... 

Case m = 1 is different from the previous one. Passage to canonical system is 
attained by means of the following transformation: 

x = U, y = W-yU, T = -VT. 

In new variables system read as follows: 

^ = -7^" (yX + y^=P2ix,y), (18) 

dy I \ fUi \ fUi V 2 / N 

— = y - — < X I — X + y 1 - Knx I— x + yl -x+UiP2{x,y) 

= y + Q2ix, y). 
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Figure 1: Phase portraits in vicinity of (0, 0) for different values of tlie parameter m 

Solving equation y + Q2{x, y) = 0, we get y = ip{x) = a2X^ + where index 02 
depends upon n Inserting function (p{x) into the RHS of the first equation we get 

P^ix, ^{x)) = A„+ia;"+^ + ... = -kUi a;"^^ + .... 

The last case we are going to analyze is that with m = |. The motivation of such a 
choice will be clear later on. In order that the analytical theory be applicable we make 
a change variables Z = \/\J . After that we get the system 

l| = _^Z2"-iH/, (19) 

^ = W{Z^~V)~ KnZ^^^-^^W^ + Z {U,- Z^) . 
Passage to the canonical variables is performed by the change of variables 

x = Z, y = W-^Z, T = -DT. 

The system resulting from this will be as follows: 

dx K _2„-l fUl 



dr 2V \V 
2 



x'^-' [^x + y] =P2{x,y), (20) 



y--Y\^yX + yj~Knx'^''-'>\^yX + yj - x'' + Ui Pi^x, y) 

= y + Q2ix, y). 

Solving first equation y + Q2{x, y) with respect to y and inserting the function '{){x) 
obtained this way into P2 we get the decomposition 

P2(x, ^{x)) = A2„x2" + ... = + - 



Using the classification given in [12] we are now able to formulate the following 
result. 
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Theorem 2. 



Ifm + n is an odd natural number then the origin of system (11) is a topological 



saddle having a pair of separatrices tangent to the vertical axis and the other pair 



tangent to the horizontal axis when m > 1 and to the line W = when m = 1. 
For even m + n the stationary point is a saddle-node with two saddle sectors lying 
in the right half-plane. Two of its three separatrices are tangent to the vertical 
axis while the remaining one is tangent to the horizontal axis in case when m > 1 
and and to the line W = when m = 1. The nodal sector lying in the left 
half-plane is unstable. 

If n is an odd number then the origin of system corresponding to / f7^) is a topo- 
logical saddle identical with that of the previous case. For even n the stationary 
point is a saddle-node identical with that of the previous case 



3. For any n>l the origin of system l[2u]) is a saddle-node geometrically identical 
with that from the previous two cases. 

The crucial fact appearing from this analysis is that the stationary points (0, 0) 
of the canonical systems f|T7|) - fl20|) depending on the values of the parameters m, n 
are either saddles or saddle-nodes with the saddle sectors placed at the right half- 
space. The return to the original coordinates does not cause the change the position 
of the saddle sectors but changes the orientation of vector fields and the angle at 
which the outgoing separatrice leaves the stationary point. The local phase portraits 
corresponding to three distinct cases are shown on figure [1] reconstructed on the basis 
if the analysis of relation between (fT7|) - (l20|) and system (fT5|) . 

Before we start to discuss the results of numerical study of system (H^]) . let us 
consider the question on when the presumably appearing homoclinic loop corresponds 
to the compacton-like solution to the source equation ([T]). To answer this question we 
are going to find the solution W{U) = ai U^^ -\- a2 U^^ + ... of equation 

KU''W^ = G = U"'(Ui-U)~Kn f/"-^ W^-(V- U) W, (21) 
dU 

(equivalent to system (fT5|) ) up to the given order and next estimate the limit of the 
vertical component of vector field 

K = li,„ GiU,W(U))^ 

Homoclinic loop will correspond to comactly-supported solution if K is finite or infinite 
and also if it tends to zero not quicker than U^, where 1/2 < /i < 1. 
Equation fl23|) can be re-written as follows: 



= [7^+^^^ (ai + a2 U^'^-^'^ + ...) -V {ai + + ...) - 



8 



The procedure of solving (123|) is pure algebraic: we collect the coefficient of different 
powers of U end equalize them to zero. The lowest power in the RHS is either f/'^i 
or If^^ or ^"+2a'i-i_ ^he number + 2 yUi — 1 cannot be less or equal to /ii because 
it involves the inequality < fii < 1 — n which is impossible for any natural n. On 
the other hand ifn + 2/ii — 1 < m then /ii becomes an "orphan" and ai should be 
nullified. 

So if^+'^t^'L-^ cannot be the lowest monomial. From this immediately appears that 
the only choice leading to the nontrivial solution is fii = m. Next important observation 
is such that the first nonzero monomial at the RHS should have as a counterpart at 
the LHS the least monomial, that is ^"+2mi-i = ^n+2m-i_ 'jj^jg turn is sufficient to 
estimate the limit fl22|) . In fact, 

G{U,W{U)) G{U,W{U)) A,U-+^^~^ + ... A, 
= = ~ — u + . . . 

A(f/) Kf/" /tf/" K 

And now we are ready to formulate the main result of this subsection. 

Proposition l.Homoclinic loop bi- asymptotic to stationary point (0, 0) of system 
( [73]) corresponds to compacton-like solution to the initial PDE if m < |. 

Presented above result delivers sufficient but not necessary condition for the homo- 
clinic loop corresponding to compacton. Basing on the local asymptotic analysis it 
is possible to show that for < m < 1 and any natural n trajectory bi-asymptotic 
to stationary point (0, 0) will also correspond to compactly-supported traveling wave. 
The proof is somewhat cumbersome and we shall not present it here. 

We would like to know, what sort of TW corresponds to the homoclinic loop in 
case when m > 1 and n > 1. Presently we are not able to answer this question 
rigorously. Nevertheless, we bring some arguments evidencing that the "tail" of TW 
in this case spreads up to — oo whereas the front sharply ends. This is so because 
the "tail" corresponds to the outgoing separatrix which is tangent to either horizontal 
axis (if m > 1) or the line W = (if m = 1). On the other hand, the incoming 
separatrix enters the origin at the angle From these observations conclusions can 
be made concerning the shape of TW. Corresponding arguments will be delivered for 
the case m = n = 1. 

Analysis of formula f|T8|) enables to state that outgoing separatrix is tangent to the 
straight line W = ^U. So the ffist coordinate of the saddle separatrix in vicinity of 
the origin satisfies the equation 

having the approximated solution f/(^) = Aexp{^^). This form of approximate solu- 
tion tells us that separatrix reaches the origin as ^ — — oo. 

Now let us consider the incoming separatrix, which is tangent to the vertical axis 
and lies at the fourth quadrant. This enables us to assume that W = —B + o (W^), 
where < a < 1 and B > 0. Approximate equation describing the ffist coordinate of 
incoming separatrix is then as follows 

^ = -BU'' + o{Un. 
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Figure 2: Phase portrait of system (fT5l) obtained for k = n = l,Ui = 3, m = ^ and 
V = 2.82 

Hence 

U{0 = [Bil-a) (^0-0]^- 
Analysis of this formula shows that trajectory reaches reaches the origin in finite "time" . 

3 Results of numerical simulation 

Numerical simulation of system (fTSjl were carried out for different values of the param- 
eters. Let us consider first the results of simulation obtained for k = n = l,Ui = 3 and 
m = |. Experiments show that at V slightly less then Vcr = f/i a stable limit cycle is 
softly created in system (ITSjl . 

Its radius grows as V decreases. As the left end of the periodic trajectory approaches 
the origin, it becomes asymmetric (fig. [2]). At V ^ 2.5397975 periodic trajectory 
destroys giving way to the homoclinic bifurcation (fig. [3]). Corresponding compacton- 
like solution to the equation ([T]) is shown on fig. |H Compacton, moving from left to 
right, is gently sloping towards the tail. It is worth noting that the compacton-like 
solution is always symmetric in case it is solution to Hamiltonian system. 

Another numerical simulations were put up for k = 1, m = 1, n = 1, f/i = 3. 
Scenario is the same as in the previous case. Limit cycle appearing a.t V < Ui grows 
as V decreases (fig. ED- Periodic trajectory is destroyed at = 2.48225, giving way 
to homoclinic trajectory (fig. E]). Corresponding TW solution to source equation ([T]) is 
shown on fig. [3 

Let us finally show patterns corresponding to stronger nonlinearity of the source 
term. We take m = 3 and the remaining parameter the same as in the previous case. 
Such a change leads to larger asymmetry of the homoclinic contour, fig. [8H9] and ap- 
pearance of wave pattern with sharp front and very long relaxing tail, reminding shock 
or detonation wave, fig. [IDl 
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Figure 5: Phase portrait of system ( JTSll obtained for K, = m = n = l, Ui = 3 and 
V = 2.775 
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Figure 6: Homoclinic bifurcation corresponding to k, = m = n = 1, Ui = 3 and 
V = 2.48225 




Figure 7: Solitary wave solution of the equation ([T]) corresponding to ip{u) 
n = l, [/i = 3 and and 1/ = 2.48225 




Figure 8: Phase portrait of system (flSll obtained for m 
V = 2.7444875 



3, K = n = 1, Ui = 3 and 
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Figure 9: Homoclinic bifurcation corresponding to m = 3, k = n = 1, Ui = 3 and 
V = 2.4449 



Figure 10: Solitary wave solution of the equation ([T]) corresponding to (p{u) = u^, 
K = n = l, f/i = 3 and and 1/ = 2.4449 
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4 Final remarks 



So we have shown that, depending on the values of the parameters convection-reaction- 
diffusion model described by equation ([T]) possesses compacton-like TW solutions or 
mixed solutions having the soliton-like tails and sharp fronts. The mere existence of 
compacton-like solutions is possible if the diffusion coefficient is a function of dependent 
variable u. In all cases considered wave patterns have been obtained as the result of 
two bifurcation. First we performed the local nonlinear analysis of dynamical system 
([5]) describing the whole set of TW solution and stated conditions making possible the 
stable limit cycle creation in proximity of critical point {Ui, 0). The other critical point, 
placed at the origin, is shown to possess simultaneously a saddle sector in the right 
half-plane, making possible the homoclinic loop appearance. Next we studied numeri- 
cally the evolutions of periodic solutions, observing how they undergo the homoclinic 
bifurcation. Analytical results obtained in subsection 12.31 enable us to state when the 
homoclinic loop obtained in numerical experiments correspond to compacton. 

The results obtained here are not completely rigorous. In particular it is rather 
unknown the precise value of wave pack velocity V corresponding to the homoclinic 
bifurcation. 

To obtain the solutions in question analytical methods based on the generalized 
symmetries [IT] and related to them anasatz-based methods [ISl [121 EOl EI] can be 
applied. Yet it is little chances to obtain them for all possible values of the parameters 
for the source equation is not completely integrable. So it would be desired to apply 
some other means, such as the rigorous computer-assisted proofs for this problem. 

In connection with the results obtained a natural question arises about the benefit 
of solutions existing merely at selected values of the parameters. So first let us men- 
tion that the velocity of wave pack, being chosen as the bifurcation parameter, is an 
"external" in some sense parameter and its value is rather connected with the portion 
of "energy" delivered at the initial moment of time, which is easily controlled, while 
the parameters characterizing system remain unchanged. 

Of prime importance is also the well known fact that invariant wave patterns such 
as kinks, solitons, compactons etc., very often play role of asymptotics, attracting in 
the long range close near-by solutions [22], [23], E] . Study of the stability and attracting 
properties of the invariant solutions to ([I]) is very non-trivial issue going far beyond 
the scope of this article. 
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